Dirac semimetal in three dimensions 
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We show that the pseudo-relativistic physics of graphene near the Fermi level can be extended to 
three dimensional (3D) materials. Unlike in phase transitions from inversion symmetric topological 
to normal insulators, we show that particular space- groups also allow 3D Dirac points as symmetry 
protected degeneracies. We provide criteria necessary to identify these groups and, as an example, 
present ab initio calculations of /3-cristobalite Bi02 which exhibits three Dirac points at the Fermi 
level. We find that /3-cristobalite Bi02 is metastable, so it can be physically realized as a 3D analog 
to graphene. 



In a Dirac semimetal, the conduction and valence 
bands contact only at discrete (Dirac) points in the Bril- 
louin zone (BZ) and disperse linearly in all directions 
around these critical points. In two dimensions, spinless 
graphene exhibits such point-like degeneracies between 
the conduction and valence bands: the low energy effec- 
tive theory at each of the critical points takes the Dirac 
form, ^^(k) = v{kx(Jx + ^y^y) where a = {cr^^, cr^, cr^} 
are the Pauli matrices and v ^ [1. The existence of 
Dirac points near the Fermi level is responsible for many 
important properties of graphene such as high electron 
mobility and conductivity. However these Dirac points 
are not robust because they can be gapped by a pertur- 
bation proportional to cr^. Spin-orbit coupling doubles 
the number of states and gaps the Dirac points [2 ; how- 
ever the splitting is very small because carbon is a light 
atom. 

In 3D, the analogous (and slightly generalized) Hamil- 
tonian is H(k) Provided det[vij] ^ 0, H{k) 

is robust against perturbations because it uses all three 
Pauli matrices. This Hamiltonian is called a Weyl Hamil- 
tonian because it describes two linearly dispersing bands 
that are degenerate at a (Weyl) point. The robustness 
of a Weyl point can be quantified by the Chern number 
of the valence band on a sphere surrounding the point, 
which takes values sgn{det[v ij]) = ±1. If a Weyl point 
occurs at some BZ momentum k, time reversal (T) sym- 
metry requires that another Weyl point occur at — k with 
equal Chern number. However, the total Chern num- 
ber associated with the entire Fermi surface must vanish. 
Hence there must exist two more Weyl points of opposite 
Chern number at k' and — k^ Inversion (I) symmetry re- 
quires that Weyl points at k and — k have opposite Chern 
number. Hence under both T and I symmetries, k = k' 
and the effective Hamiltonian involves four linearly dis- 
persing bands around k. Such a Hamiltonian is called a 
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Dirac Hamiltonian, and it is not robust against pertur- 
bations because there are additional 4x4 Dirac matrices 
that can be used to open a gap at the Dirac point. 

The Fermi surface of a Dirac semimetal consists en- 
tirely of such point-like (Dirac) degeneracies. 3D Dirac 
semimetals are predicted to exist at the phase transi- 
tion between a topological and a normal insulator when 
I-symmetry is preserved [31 [4] (Ref. [5] demonstrates 
such a Dirac point degenerate with massive bands.) If 
either I or T-symmetry is broken at the transition, a 
Dirac point separates into Weyl points and one ob- 
tains a Weyl semimetal (Fig. |l(c)"] ). The topological na- 
ture of Weyl points gives rise to interesting properties 
such as Fermi- arc surface states [6] and pressure induced 
anomalous Hall effect 0. Recent proposals to design 
a Weyl semimetal have been predicated upon the exis- 
tence of a parent Dirac semimetal which splits into a Weyl 
semimetal by breaking I [8 or T-symmetry [9]. Ref. [10] 
demonstrates the existence of bulk chiral fermions due to 
crystal symmetry in single space-groups. 

Dirac points that arise in a topological phase transition 
described above are accidental degeneracies. In general, 
two Weyl points with opposite Chern numbers annihi- 
late each other unless their degeneracy is protected by 
additional space-group symmetry. Therefore we ask if a 
Dirac point can arise as a result of a crystallographic sym- 
metry. Indeed certain double space-groups allow Dirac 
points at high symmetry points on the boundary of the 
BZ. As an example we present ab initio calculations of /3- 
cristobalite Bi02 (Fig. |2(b)[ ) which exhibits Dirac points 
at three symmetry related X points on the boundary of 
the FCC BZ (Figs. [T(b)] and [3(c)| . This system realizes 
a Dirac degeneracy first encountered in a tight-binding 
model of s-states in diamond in Ref. jTT] (the Fu-Kane- 
Mele model). In the absence of T-symmetry, two Weyl 
points with equal Chern numbers can be degenerate due 
to a point group symmetry as shown in Ref. [12 . 

A 3D double space-group must satisfy the following 
criteria to allow a Dirac point. It must admit four di- 
mensional irreducible representations (FDIRs) at some 
point k in the BZ such that the four bands degenerate 
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at k disperse linearly in all directions around k and the 
two valence bands carry zero total Chern number. If 
the little group Gk at k contains a three-fold or a six- 
fold rotation symmetry and the valence and conduction 
bands around k are non-degenerate, the Chern number 
of the FDIR is guaranteed to be non-zero. This rules out 
symmorphic space-groups with FDIRs because they con- 
tain three-fold rotations. This also rules out interior BZ 
momenta because non-symmorphic little groups without 
three-fold rotations exhibit FDIRs only on the boundary 
of the BZ [13 . To guarantee linear dispersion of bands 
around k, the symmetric kronecker product [R\^ x R\^] 
of the FDIR with itself must contain the vector repre- 
sentation of Gk [14 . Finally, away from k, the FDIR 
must split so that the valence and conduction bands are 
non-degenerate everywhere except at k (Fig. [4| . 

We apply the above criteria to two important space- 
groups. The space-group of diamond (227, Fd3m), which 
is also the symmetry group of /3-Cristobalite Bi02, ex- 
hibits FDIRs Rr at F and Rx at X. Gr contains three- 
fold rotation symmetry and [Rr x Rr] does not contain 
the vector representation of Gr- Therefore the F point 
in a diamond lattice cannot host a Dirac point. Rx is 
a projective representation of Gx which does not have 
any three-fold rotations because all the point group op- 
erations in Gx are those of the group D^h- [Rx x Rx] 
contains the vector representation of Gx- Finally Rx 
splits into either two doublets or four singlets away from 
X (Figs.[4];a) andgb)). Therefore the X point m space- 
group 227 is a candidate to host a Dirac semimetal if its 
FDIR can be elevated to the Fermi level. Indeed we show 
that /3-Cristobalite Bi02 exhibits such a Dirac point at 
X, Fig. |3(c)| The Dirac point at X in the FK M mo del 
is also spanned by states belonging to Rx (Fig. |3(d) ). 

The zincblende lattice (space-group 216, F43m) has an 
FDIR at F and the little group Gp has a three-fold ro- 
tation symmetry. [R[. x i?p] contains the vector represen- 
tation of Gp. Mirror symmetry in Gp requires i?p to split 
into a two-fold degenerate representation and two non- 
degenerate representations along the (111) axis, which is 
also the symmetry axis for the three-fold rotation. Time 
reversal symmetry requires that the two-fold degenerate 
band remain flat along the (111) axis. Fig. ^d). Thus 
the lowest band carries Chern number 0, while the two 
fiat bands carry 1 and -1. Therefore the dispersion of i?p 
is not Dirac-like along (111). 

In HgTe, which takes the zincblende lattice, the degen- 
erate valence and conduction states at F span i^p and 
constitute the entire Fermi surface. It is known that in 
HgTe the valence and conduction bands disperse linearly 
in two directions around F and quadratically in a third 
(Fig. [4]^d) and Ref. [TF). One might ask if a perturba- 
tion might turn HgTe into a Dirac semimetal. However 
the zincblende lattice does not satisfy the criteria for 3D 
Dirac points as outlined above, so HgTe cannot host a 
Dirac semimetal. (a) F is an interior point of the BZ and 



the little group at F contains a three-fold rotation, (b) 
Mirror symmetry requires two bands to be degenerate 
along the (111) axis but since the Chern number must 
vanish, the degenerate bands must be fiat and consist of 
a conduction and a valence band. This is why we see 
quadratic dispersion along the (111) axis, (c) Breaking 
mirror symmetry splits the degenerate flat band but then 
the Fermi surface develops other non-Dirac like pock- 
ets to compensate for the non-zero Chern number, (d) 
Breaking three-fold rotation symmetry splits the degen- 
eracy at F entirely and the material becomes a topologi- 
cal insulator [16]. 

We briefly discuss the theory behind the above crite- 
ria. We are interested in FDIRs of double space-groups 
at points k such that the valence and conduction bands 
are distinct in a small region around k and carry zero 
total Chern number. The Chern number of a degenerate 
representation can be determined up to an integer by the 
rotation eigenvalues of the valence bands. Electron states 
spanning an FDIR are equivalent to a quadruplet 
which exhibit eigenvalues g±*37r/n g±27r/n a 27r/n ro- 
tation symmetry. Rotation eigenvalues of states at time 
reversed momenta about the degenerate point are com- 
plex conjugates. Therefore the FDIR will carry Chern 
numbers ±1 mod n for one valence band and ±3 mod n 
for the other with total Chern number ±4 mod n or ±2 
mod n for the FDIR. This is zero only for n = 1, 2, 4. If 
the conduction and valence bands are distinct in a small 
region around k, the Chern number of the FDIR will be 
non-zero if the little group Gk contains a 27r/3 or 27r/6 
rotation symmetry. In HgTe however, the little group at 
r contains a three-fold rotation symmetry but the FDIR 
at F has zero Chern number because one of the valence 
bands is degenerate with one of the conduction bands 
along the (111) axis. 

Non-symmorphic space-groups contain point group op- 
erations coupled with non-primitive lattice translations. 
For example, inversion interchanges the FCC sub-lattices 
in the diamond space-group. Representations of non- 
symmorphic space-groups at momenta inside the BZ mo- 
menta are obtained from regular representations, while 
those at the surface BZ momenta are obtained from 
projective representations of the associated crystal point 
group. The factor system of the projective representa- 
tion is chosen to implement the required non-primitive 
translation corresponding to the non-symmorphic point 
group operation [13]. A theorem by Schur guarantees 
that projective representations of a group can be ob- 
tained by restricting to the group elements the regular 
representations of a larger group called the central exten- 
sion group \TT . The central extension of a group is ob- 
tained by taking its product with another finite Abelian 
group. The important point to emphasize is that repre- 
sentations of non-symmorphic space-groups are obtained 
from representations of central extensions of the 32 point 
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groups. Central extension groups exhibit FDIRs even 
without three- fold rotations in the original point group. 
This is precisely why Dirac points can exist in 3D as 
symmetry allowed degeneracies. 

To realize a Dirac-like dispersion in the vicinity of an 
FDIR, some of the matrix elements {ipi\p\tljj)^ where 
span the FDIR, must be non-zero. This is guaranteed if 
the symmetric kronecker product of the FDIR with itself 
contains the vector representation of the central exten- 
sion group to which the FDIR belongs J4]. We restrict 
to the symmetric part of the kronecker product because 
matrix elements {ipi\p\tljj) correspond to level transitions 
between states spanning the same representation Fi- 
nally, the allowed representations in the vicinity of an 
FDIR should be such that each band disperses with non- 
zero slope in all directions. This is possible only if the 
valence band is distinct from the conduction band ev- 
erywhere except at the Dirac point. Fig. [4] illustrates the 
various possible ways in which an FDIR can split linearly. 

Although crystallographic symmetries determine 
whether 3D Dirac points can exist, physical and chem- 
ical considerations dictate whether they arise at the 
Fermi level without additional non-Dirac like pockets in 
the Fermi surface. In the FKM model, the Dirac point 
at X appears at the Fermi energy. However, in known 
materials on a diamond lattice 5-states appear below the 
Fermi energy. In realistic systems, additional orbitals 
hybridize with these s-states and bands cross the Fermi 
level at other points besides X. The problem is espe- 
cially severe in space-group 227: without spin, the line 
V from X to W is two-fold degenerate. With spin-orbit 
coupling, this line splits weakly for lighter atoms so 
the bands dispersing along this line can hybridize and 
introduce additional Fermi surface. Forcing species with 
valence states on the diamond lattice would fail to 
realize the FKM model. Indeed, ab initio calculations 
with group I elements and gold show that the splitting 
along V is insufficient to overcome this dispersion. In 
some cases, additional bands crossed the Fermi level. We 
performed ab initio calculations using the plane wave 
density functional theory package quantum espresso [18^ , 
and designed non-local pseudopotentials [191 EH] with 
spin-orbit interaction generated by OPIUM. 

We consider derivatives of the diamond lattice that 
remain in space-group 227. We place additional atoms 
in the lattice such that the configuration of added species 
allows its valence orbitals to either belong to the FDIR 
of interest, or appear away from the Fermi energy of the 
final structure. If the new species can split the nearby p 
states of the existing atoms away from the s levels, band 
crossing at the Fermi level can be avoided. 

One such structure is /3-cristobalite Si02 (Fig. |2(b)] ), 
which consists of silicon atoms on a diamond lattice with 
oxygen atoms placed midway along each silicon-silicon 
bond [IT]- Oxygen atoms have two consequences: part 
of the O p-shell strongly hybridizes with the Si p-states. 



moving them away from the Si 5-states, while the re- 
maining O p-states span the same representation as the 
Si 5-states. A Dirac point can be realized by Si s—0 p 
bonding/anti-bonding set of states. Fig. 2(a) shows that 
the Si 5—0 p bands are present and take a configuration 
similar to the valence and conduction bands in the FKM 
model, but appear well below the Fermi energy. Addi- 
tionally, the bands are nearly degenerate along the line 
V from X to W due to weak spin-orbit coupling. 

Heavier atoms substituting Si both widen this gap and 
bring the FDIR of interest at X to the Fermi level. Fig. [3] 
shows the band structures of compounds /3-cristobalite 
XO2 where X = As/Sb/Bi. The change in chemical iden- 
tity promotes the X s—0 p four- fold degeneracy at X to 
the Fermi level, and stronger spin-orbit coupling widens 
the gap along V. Bi02 bears striking similarity to the 
FKM model, with linearly dispersing bands in a large 
energy range around a Dirac point at the Fermi level. 
Our calculations show that the phonon frequencies for /3- 
cristobalite Bi02 at F are positive, so it is a metastable 
structure. Further calculations reveal that it becomes un- 
stable under uniform compression exceeding 2GPa, which 
represents a stability barrier of approximately 0.025eV 
per atom. On this basis, the possibility of synthesis ap- 
pears promising. However, Bi204 is also likely to take 
the cervantite structure (after Sb204, which has similar 
stoichiometry [21]) which is 0.5 eV per atom lower in en- 
ergy as compared to /3-cristobalite and 60% smaller in 
volume. Therefore we conclude that /3-cristobalite Bi02 
would be metastable if synthesized, although preventing 
it from directly forming the cervantite structure would be 
challenging. Nonetheless we have provided an existence 
proof of a Dirac semimetal in /3-cristobalite Bi02 due to 
real atomic potentials at the DFT level. 
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FIG. 1: 3D Dirac semimetal in /3-cristobalite Bi02. (a) Bril- 
louin zone (BZ) of the FCC lattice. The plane highlighted 
in gray joins the three symmetry related X points. Other 
high symmetry points are also indicated. |(b)| Conduction and 
valence bands of /3-cristobalite Bi02 are plotted as functions 
of momentum on the plane highlighted in gray on the left. 
Each band is two-fold degenerate due to inversion symme- 
try. Dirac points appear at the center of the three zone faces 
of the BZ. |(c)| Dirac, Weyl and insulating phases in the dia- 
mond lattice. (1) The states at the Dirac point at X span a 
four dimensional projective representation of the little group 
at X which contains a four-fold rotation accompanied by a 
sub-lattice exchange operation. (2) Four Weyl points on the 
zone face due to a small inversion breaking perturbation. The 
Chern number of each Weyl point is indicated. (3) Two Weyl 
points appear on the line from XioW for a T-breaking Zee- 
man field B oriented along that direction. B oriented along 
other directions gaps all the Dirac points by breaking enough 
rotational symmetry that no two-dimensional representations 
are allowed. (4) Gapped phase obtained by breaking the four- 
fold rotation symmetry or by applying a magnetic field in any 
direction except along x, ^, or z. The insulating phase can be 
a normal, strong or a weak topological insulator [fTj. 



(a) (b) 
FIG. 2: (a) Band structure of /3-cristobalite Si02. Energy 



bands are plotted relative to the Fermi level. Each band is 
two-fold degenerate due to inversion symmetry. The (high- 
lighted) FDIR at —4.5 eV is split into two linearly dispersing 
bands between X and F while the two degenerate bands along 
X and W are weakly split. This FDIR is buried deep be- 



low the Fermi level, (b) The /3-cristobalite structure of Si02 
(Bi02). Silicon (bismuth) atoms (light gray) are arranged 
on a diamond lattice, with oxygen atoms (dark gray) sitting 
midway between pairs of silicon (bismuth). 
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FIG. 3: Band structures of[(aJlAs 02,|(b)1 Sb02, and[(^BiO 



in the /3-cristobalite structure, and | (d) | instates on a diamond 
lattice in the tight-binding model of Ref. 11 . Energy bands 
are plotted relative to the Fermi level. Each band is two-fold 
degenerate due to inversion symmetry. Insets: with increas- 
ing atomic number of the cation, spin-orbit coupling widens 
the gap along the line V from X to W. In Bi02 and Sb02, the 
dispersion around the X point is linear in all directions indi- 
cating the existence of Dirac points at X. Bi02 and Sb02 are 
Dirac semimetals because their Fermi surface consists entirely 
of Dirac points. 
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FIG. 4: Linear splitting of four-fold degenerate irreducible 
representations (FDIRs). If the symmetric kronecker prod- 
uct of an FDIR with itself contains the vector representation 
of the group to which the FDIR belongs, it will split in one 
of the four possible ways displayed above, (a) The FDIR 
splits into two two- fold degenerate bands. This situation is 
realized at the X point of the FCC Brillouin zone in a dia- 
mond lattice, (b) The FDIR splits into four non-degenerate 
bands. This situation arises at the F point in zincblende if 
mirror symmetry is broken (although the FDIR in zincblende 
develops a non-zero Chern number due to three-fold rotation 
symmetry at F). (c) The FDIR splits into two non-degenerate 
and one two-fold degenerate band with linear dispersion, (d) 
The splitting of the FDIR at F in zincblende. The two-fold 
degenerate band is constrained to be flat implying quadratic 
dispersion along that direction. The Chern number of this 
representation is zero inspite of a three-fold rotation symme- 
try because the conduction and valence bands are degenerate 
away from F. 



